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Abstract 

Using ab initio calculations based on density-functional theory we have performed a theoretical 
investigation of substitutional boron and nitrogen impurities in silicon carbide (SiC) nanoribbons. 
We have considered hydrogen terminated SiC ribbons with zigzag and armchair edges. In both 
systems we verify that the boron and nitrogen atoms energetically prefer to be localized at the 
edges of the nanoribbons. However, while boron preferentially substitutes a silicon atom, nitrogen 
prefers to occupy a carbon site. In addition, our electronic-structure calculations indicate that 
(i) substitutional boron and nitrogen impurities do not affect the semiconducting character of the 
armchair SiC nanoribbons, and (ii) the half-metallic behavior of the zigzag nanoribbons is main- 
tained in the presence of substitutional boron impurities. In contrast, nitrogen atoms occupying 
edge carbon sites transform half-metallic zigzag nanoribbons into metallic systems. 
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I. INTRODUCTION 



Silicon carbide (SiC) nanostructures have attracted much interest in the last years be- 



cause they combine the excellent properties of SiC bulk[l|-|3| (wide band gap, high thermal 
conductivity, high breakdown electric field, thermal stability, high electronic mobility, and 
others) with quantum-size effects. This combination makes SiC nanostructures interesting 
materials for nanotechnology applications and nanoscale engineering, 

Recently, increasing attention has been paid to the planar structures of SiC. Although 
silicon prefers sp^ instead of sp^ hybridization, SiC nanoribons (SiCNRs) have been suc- 
cessfully synthesized jsl and first-principles calculations have shown that SiC has a stable 
two-dimensional honeycomb structure, similar to graphene. {g-S] In addition, recent ab ini- 
tio studies revealed that the monolayer SiC sheet is semiconducting with a band gap of about 
2.55 eV,^-8| while armchair SiCNRs are nonmagnetic semiconductors (with a wide band 

nn 

gap for all widths), |9|,ll0l and zigzag SiCNRs narrower 4 nm present half-metallic behavior 
without applied external electric field or chemical modifications, [q 



11| This feature makes 



narrow zigzag SiCNRs promising candidates for spintronic applications. 

The great potential of the SiC nanoribbons for future applications in spintronics stimu- 
lates detailed studies of the effects of impurities and defects in their electronic properties. 
In particular, SiC nanoribbons doped with boron and nitrogen atoms is an important issue 
to be investigated, since B and N are common contaminants in SiC bulk, acting as p- and n- 
type dopants, respectively. Moreover, it is well known that substitutional boron and nitrogen 
impurities in zigzag graphene nanoribbons break the spin degeneracy of the transmittance 
channels, and transform metallic nanoribbons into semiconducting systems. 12|-|l5| 



In this work we have performed a theoretical investigation, using ab initio calculations 
based on the density-functional theory, of substitutional boron and nitrogen impurities in 
armchair and zigzag SiC nanoribbons. Our results indicate that the boron and nitrogen 
atoms energetically prefer to be localized at the edges of both nanoribbons. However, while 
boron preferentially substitutes a silicon atom, nitrogen prefers to occupy a carbon site. In 
addition, we verify that substitutional B and N impurities do not affect the semiconducting 
character of the armchair SiCNRs, although the spin degeneracy of the band structure 
is broken. On the other hand, we observe different effects of these impurities on the half- 
metallic behavior of the zigzag SiCNRs: the half-met allicity of the nanoribbons is maintained 



in the presence of substitutional B impurities, whereas zigzag SiCNRs become metalhc when 
doped with N atoms. 



II. METHOD OF CALCULATION 



The ab initio calculations presented in this work were carried out within the frame- 
work of the density functional theory (DFT), 16 1 using the SIESTA code.[l7| We employed 
the generalized grad ient approximation as implemented by Perdew, Burke and Ernzer- 
hof (GGA-PBE) 181] for the exchange-correlation functional, and norm-conserving fully- 
separable pseudopotentials 19j to treat the electron-ion interactions. The Kohn-Sham or- 



20| using 



bitals were expanded in a linear combination of numerical pseudoatomic orbitals 
a split- valence double-zeta basis set with polarization functions (DZP). 21] 

We have considered hydrogen-passivated zigzag and armchair SiCNRs with widths W = 
4 (4-ZSiCNR) andW = 7 (7-ASiCNR),|22| respectively. The 4-ZSiCNR and 7-ASiCNR 
structures were modeled within the supercell approach with 40 and 54 atoms (including 
hydrogen atoms) in the unit cell, respectively. The structural models of 4-ZSiCNR and 
7-ASiCNR, indicating the substitutional sites investigated, are shown in Fig. [H A vacuum 
region of about 20 A along the nonperiodic directions was employed to avoid interactions 
between two neighboring ribbons. All the atomic positions were relaxed during the geometry 
optimization until the Hellman-Feynman forces were below 0.01 eV/A. The Brillouin zone 
was sampled using the Monkhorst-Pack scheme 231] with a (11 x 1 x 1) grid for the total- 
energy calculations and a (120 x 1 x 1) mesh for the electronic-structure calculations. The 
convergence of our total-energy results with respect to the number of k-points was verified. 

In each system the energetic stability of the impurities was determined by comparing 
their formation energies 2d, 125] 

n = ET-Y.N^^i^^ (1) 



where is the total energy of the structure, fii is the chemical potential of the atomic 
specie i, and Ni is the number of atoms i in the structure. The chemical potentials of Si 
and C are restricted to the ranges 



- AH{S\ 



(2) 
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FIG. 1. (Color online) Structural models of (a) 4-ZSiCNR and (b) 7-ASiCNR indicating the 
substitutional sites investigated. Odd (even) numbers correspond to Si (C) sites. 



and 



- (SiC) + //^"^'^ < /ic < /ig 



.bulk 



(3) 



where Ai^f(SiC) is the SiC formation heat. 

The (relative) formation energies presented in this work were calculated assuming Si/C 
stoichiometric condition 



/^c,si = /iS-^A^(SiC). 
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(4) 



III. RESULTS AND DISCUSSIONS 



Initially we examined the equilibrium geometries and the electronic structures of pristine 
4-ZSiCNR and 7-ASiCNR. We found optimized Si-C bond lengths from 1.77 to 1.84 A in 
4-ZSiCNR, and from 1.75 to 1.81 A in 7-ASiCNR. In both systems, the Si-H and C-H bond 



lengths were about 1.52 and 1.11 A, respectively. These values are in goo d a, 



26 



greement with 



27| By analysing 



previous calculations for SiC nanoribbonsJlO'] sheets[6j and nanotubes. 
the electronic structure in Fig. [2] we observed that pristine 4-ZSiCNR [Fig. |2] (a)] presents a 
half-metallic behavior, with a semiconducting spin-up channel (band gap of 0.22 eV at the F 
point) and a metallic spin-down channel. On the other hand, pristine 7-ASiCNR [Fig. |2] (b)] 
is a semiconductor with a direct band gap of 2.22 eV (at the F point) for both spin-up and 
spin-down channels. 



FIG. 2. (Color online) Electronic band structures of pristine (a) 4-ZSiCNR and (b) 7-ASiCNR. The 
band structure of pristine 4-ZSiCNR is depicted in two different ranges: (al) \E — Ef\ < 3.0 eV 
and (a2) \E — Ef\ < 0.2 eV. Spin-up and spin-down bands are indicated by solid red and dashed 
blue lines, respectively. Dotted line indicates the position of the Fermi level. 

Considering one B or N atom per supercell we investigated 8 substitutional sites (sites 
labeled from 1 to 8 in Fig. [1]) in both 4-ZSiCNR and 7-ASiCNR. In B-doped nanoribbons 
(B/4-ZSiCNR and B/7-ASiCNR) the formation energy of each B\ configuration was com- 
pared with that of Bg^ configuration 



Here i G {!,..., 8}, A G {Si, C} and Bg^ denotes the configuration in which B substitutes 
Si at the site 1. In the N/4-ZSiCNR and N/7-ASiCNR systems, we compared each 
configuration with configuration using 



Our results for the relative formation energies AQb and A^n are summarized in Table [B 
In this case, positive values of A^b (A^n) indicate that Bg^ (N|.) is energetically more 
favorable than B\ (N^). 

As can be seen in Table I, Bg; and are the most favorable configurations in both 
4-ZSiCNR and 7-ASiCNR. This indicates that the B and N atoms energetically prefer to 
be localized at the egde of the nanoribbons. However, while B preferentially occupies a Si 
site, N prefers to substitute a C atom. Among the Be configurations, B^ and B^ are the 



A^B = n{B\) -n{Bl^). 



(5) 



A^N = fi(NX) -fi(N^). 



(6) 
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most stable in 4-ZSiCNR and 7-ASiCNR, respectively. However, at the Si/C stoichiometric 
condition, (B^) is energetically less favorable than Bg; by 1.16 (0.89) eV in 4-ZSiCNR 
(7-ASiCNR). Extending the calculations of the formation energies Q(Bq), ^(Bq) and f2(Bgi) 
for different values of fic and fisi, we verify that ^^(B^) > i7(Bgj) and f^(Bc)> f2(Bgj) within 
the limits given by the expressions ([2]) and (E]) with Ai7(SiC) = 0.72 eV. |28| These results 
indicate that Be is not expected to occur in both armchair and zigzag SiCNRs. Performing 
the same analysis for N impurities, we observe that Ngj is the most likely configuration 
among the Ngi configurations in both 4-ZSiCNR and 7-ASiCNR, but Ngj is less favorable 
than within the allowed ranges for iiq and fisi (expressions ([2]) and ([3])). Thus, it is not 
expected that N atoms occupy Si sites in SiC nanoribbons. 

The energetic preference of the B (N) atoms for Si (C) sites at the stoichiometric condition 



has been also observed in SiC nanotubes 



29| and nanowires. 30| In these nanostructures. 



like in SiC nanoribbons, the formation of Ngi is not expected to occur. Nevertheless, in 
SiC nanotubes and nanowires Be becomes more likely than Bsi at the limit ^si = l^sf^^ 
= ^buik _ A/f(SiC).|29|, bo] 

We next examined the equilibrium geometries of the most likely configurations (Bg; and 
Nq) in 4-ZSiCNR and 7-ZSiCNR. The variations in the Si-C bond lengths with respect 
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FIG. 3. (Color online) Equilibrium geometries of (a) B^i/4-ZSiCNR, (b) N^/4-ZSiCNR, (c) 6^7- 
ASiCNR, and (d) Ny7-ASiCNR. The atomic distances are in A. 

to the pristine structures were less than 0.03 A in the By4-ZSiCNR, Ny4-ZSiCNR and 
B^i/7-ASiCNR systems, and less than 0.05 A in Ny7-ASiCNR. As depicted in Fig. El we 
found (i) B-C and B-H bond lengths of 1.57 and 1.24 A (1.54 and 1.24 A), respectively, 
in the By4-ZSiCNR (B^i/7-ASiCNR) system, and (ii) N-Si and N-H bond lengths of 1.75 
and 1.03 A, respectively, in Ny4-ZSiCNR. In the Ny7-ASiCNR structure, besides the N-Si 
bond lengths of 1.72 and 1.80 A, we also observed that the Sii atom [see Fig. (U^b)] moves 
upward (along the z direction) by 0.28 A with respect to the pristine 7-ASiCNR, while the 
hydrogen atom bonded to Sii is displaced by 0.77 A in the opposite direction. A similar 
behavior is found in NQ/7-ASiCNR (the second more stable configuration for N/7-ASiCNR), 
where the Sii and H atoms are displaced by 0.25 and 0.87 A, respectively. By constraining 
the relaxation of the NQ/7-ASiCNR and N^/7-ASiCNR systems along the z direction (the 
atoms are free to relax in the x and y directions), we verify that their total energies are 
increased by 0.60 and 0.62 eV, respectively, in comparison with the fully relaxed structures. 
Even so, remains the most favorable configuration for N/7-ASiCNR. 

Figures m and [5] display the electronic band structures and the projected density of states 
(PDOS) of By4-ZSiCNR, Ny4-ZSiCNR, B^ 7-ASiCNR, and 7-ASiCNR. In Fig. g^a) 
we observe that the Bg-/4-ZSiCNR system presents half- metallic characteristics: the spin- up 
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FIG. 4. (Color online) Electronic band structures of (a) B^i/4-ZSiCNR, (b) Ny4-ZSiCNR, (c) 
Bg;/7-ASiCNR, and (d) Nq/T-ASICNR. Spin-up and spin-down bands are indicated by solid red 
and dashed blue lines, respectively. Dotted line indicates the position of the Fermi level. 

channel is semiconducting, with a direct band gap of 1.30 eV at the F point, and the spin- 
down channel is metallic, with an energy band crossing the Fermi level (Ep) near to the X 
point. This result indicates that one substitutional B impurity per supercell does not affect 
the half-metallic behavior of 4-ZSiCNR. The PDOS of the By4-ZSiCNR system [Fig. [5](a)] 
reveals that the partially occupied spin-down electronic state, as well as the highest-occupied 
spin-up state, is mostly composed by C 2p and Si 3d orbitals, without contributions from 
Bgj. The B 2p orbitals are localized below {Ep — 0.5) eV and above {Ep + 0.9) eV. 

As depicted in Figs. |l](b) and[5]^b), the N|;/4-ZSiCNR system is metallic, with degenerate 
spin-up and spin-down channels. The partially occupied state comes mainly from the Si 3p 
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FIG. 5. (Color online) Projected density of states (PDOS) of (a) B^i/4-ZSiCNR, (b) Ny4-ZSiCNR, 
(c) B^i/7-ASiCNR, and (d) Ny7-ASiCNR. Solid thick red lines indicate the PDOS of the B or 
N atoms. Dashed blue and solid (thin) black lines represent the PDOS of the Si and C atoms, 
respectively. The position of the Fermi level is indicated by dotted lines. Positive and negative 
values of PDOS correspond to spin-up and spin-down contributions, respectively. Left vertical 
scale: PDOS of Si and C atoms. Right vertical scale: PDOS of B or N atoms. 

orbitals localized along the zigzag edge opposite to N|.. The occupied state situated between 
1.2 and 1.4 eV below the Fermi level is mostly composed by Si 3p, Si 3d, C 2p and N 2p 
orbitals. 

Figures SJ^c) and ID^d) show that the semiconducting character of pristine 7-ASiCNR is 
maintained in the presence of substitutional B and N impurities, but the spin degeneracy 
of the band structure is broken. The Bgj/T-ASiCNR system [Fig. IH^c)] presents a direct 
band gap of 2.17 eV (at F point) for the spin-up channel and an indirect band gap of 
0.39 eV for the spin-down channel. We found an unoccupied spin-down state lying between 



0.1 and 0.3 eV above the Fermi level. This state consists mainly of B 2p, C 2p, Si 3p and 
Si 3d orbitals [Fig. [5t^c)]. In N^/T-ASiCNR [Fig. Ht^d)], both spin-up and spin-down channels 
exhibit semiconducting characteristics, with direct band gaps of 1.34 eV (at the X point) and 
2.15 eV (at the F point), respectively. The highest-occupied spin-up state is very localized 
and lies at around 0.6 eV below the Fermi level. The PDOS of Ny7-ASiCNR [Fig. EJ^d)] 
indicates that this state is mainly composed by Si 3p, Si 3d and N 2p orbitals, and by the 2p 
orbitals of the inner C atoms. We verified that the C 2p orbitals localized along the edges 
do not contribute to this electronic state. 



IV. CONCLUSIONS 



In summary, we performed a theoretical investigation, using ab initio calculations, of 
substitutional B and N impurities in SiC nanoribbons. In both ASiCNR and ZSiCNR 
we found an energetic preference for B and N atoms occupying edge sites. However, we 
verified that B preferentially substitutes a Si atom, whereas N prefers to occupy a C site. 
We also observed that the formation of Be and Nsi is not expected to occur in SiCNRs. 
In addition, our electronic-structure calculations revealed that B and N impurities have 
different effects on the electronic character of SiCNRs. We observed that (i) substitutional 
B and N impurities do not affect the semiconducting behavior of ASiCNRs, and (ii) the half- 
metallicity of ZSiCNRs is maintained in the presence of substitutional B atoms. In contrast, 
we verified that half-metallic ZSiCNRs become metallic when doped with N atoms. These 
results suggest that the electronic properties of SiCNRs can be controlled by B- and N-doping 
processes. 
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